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We construct approximate initial data for non-spinning black hole binary systems by asymptoti- 
cally matching the 4-metrics of two tidally perturbed Schwarzschild solutions in isotropic coordinates 
to a resummed post-Newtonian 4-metric in ADMTT coordinates. The specific matching procedure 
used here closely follows the calculation in 0, and is performed in the so called buffer zone where 
both the post-Newtonian and the perturbed Schwarzschild approximations hold. The result is 
that both metrics agree in the buffer zone, up to the errors in the approximations. However, since 
isotropic coordinates are very similar to ADMTT coordinates, matching yields better results than in 
the previous calculation where harmonic coordinates were used for the post-Newtonian 4-metric. 
In particular, not only does matching improve in the buffer zone, but due to the similarity between 
ADMTT and isotropic coordinates the two metrics are also close to each other near the black hole 
horizons. With the help of a transition function we also obtain a global smooth 4-metric which 
has errors on the order of the error introduced by the more accurate of the two approximations we 
match. This global smoothed out 4-metric is obtained in ADMTT coordinates which are not horizon 
penetrating. In addition, we construct a further coordinate transformation that takes the 4-metric 
from global ADMTT coordinates to new coordinates which are similar to Kerr-Schild coordinates 
near each black hole, but which remain ADMTT further away from the black holes. These new 
coordinates are horizon penetrating and lead, for example, to a lapse which is everywhere positive 
on the t = slice. Such coordinates may be more useful in numerical simulations. 

PACS numbers: 04.20.Ex, 04.25. Dm, 04.25. Nx, 04.30.Db, 95.30.Sf 



I. INTRODUCTION 

The modeling of binary black hole systems is essential 
for the detection of gravitational waves both by space 
and earth-based interferometers [3- Q3- C3- [3- II3| - Often 
this detection relies on matched filtering, which consists 
of comparing accurate waveform templates to the signal. 
In the strong field regime, these templates can be mod- 
eled accurately enough only by numerical simulations , 
where one solves the Einstein equations subject to ini- 
tial conditions that determine the physical character of 
the system. The astrophysical accuracy of the template 
generated via a numerical evolution will depend on two 
main factors: the numerical error introduced by the evo- 
lution and gravitational wave extraction codes; and the 
astrophysical accuracy of the initial data set. In this pa- 
per we study one particular approach to solve the initial 
data problem. This approach was developed in Ref. 1] 
(henceforth paper 1) and is based on matching several 
approximate solutions of the Einstein equations. 

There have been numerous efforts to find astrophysi- 
cally accurate initial data for binar y sy ste ms l8 |. [ gl, |l ol. f ill 
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which usually consist of the 3-metric and extrinsic curva- 
ture. Most methods rely on performing decompositions 
of the data to separate the so called physical degrees of 
freedom from those constrained by the Einstein equa- 
tions and those associated with diffeomorphisms. These 



physical degrees of freedom carry the gravitational wave 
content of the data, and thus should be chosen to rep- 
resent the astrophysical scenario that is being described. 
However, the exact solution to the Einstein equations for 
the 2-body scenario remains elusive. Therefore, the grav- 
itational wave content of the data is not known exactly. 

One commonly used cure to this problem is to give 
an ansatz for these physical degrees of freedom, which 
usually does not satisfy the constraints. This ansatz is 
then projected onto the constraint-satisfying hypersur- 
face by solving the constraint equations for the quan- 
tities introduced by the decomposition. The projected 
data will now satisfy the constraints, but it is unclear 
whether such data are astrophysically accurate. First, 
some of the assumptions that are used in the construc- 
tion of the initial ansatz are known to be physically inac- 
curate. For example, the spatial metric is often assumed 
to be conformally flat, which is known to be incorrect at 
2 post-Newtonian (PN) order [2^, and which also docs 
not contain realistic tidal deformations near the black 
holes. Furthermore, after projecting onto the constraint 
hypersurface, the data will be physically different from 
the initial ansatz 0, ITiL |29^ | . It is thus unclear exactly 
what physical scenario these projected data represent and 
whether or not they will be realistic. 

Recent efforts have concentrated on using physical ap- 
proximations to construct initial data sets. One such 
effort can be called the quasi-equilibrium approach 0, 
I20L [2(| where one assumes that the quantities describ- 
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ing the initial data evolve on a time scale much longer 
than the orbital timescale, when appropriate coordinates 
are used. There are also approaches built entirely from 
second pos t-Newtonian approximations [TtI l3C| . Even 
though [T3 also develops a complete method to project 
onto the constraint hypersurface, both |Tt| and |3fT ] are 
not guaranteed to be realistic close to the black holes, 
since post-Newtonian theory is in principle not valid close 
to black holes. Another recent effort, which will be the 
subject of this work, involves several analytic approxi- 
mations that are both physically accurate and close to 
the constraint hypersurface far and close to the holes 
(see paper 1). Such initial data coming from analytic 
approximations can then be evolved without solving the 
constraints, provided that constraint violations are ev- 
erywhere smaller than numerical error. A perhaps more 
appealing alternative is to project these data onto the 
constraint hypersurface. If these data are close enough 
to the hypersurface, then any sensible projection algo- 
rithm should produce constraint-satisfying data that are 
close to the original approximate solution. 

Regardless of how such data are implemented, there 
are reasons to believe that such an approach will pro- 
duce astrophysically realistic initial data. First, since 
the data are built from physical approximations, there 
are no assumptions that are unrealistic. The only limit 
to the accuracy is given by the order to which the ap- 
proximations are taken. Moreover, the analytic control 
provided by these physical approximations allows for the 
tracking of errors in the physics, thus providing a means 
to measure the distance to the exact initial data set. Note 
however, that our work described below will only be ac- 
curate up to leading order in both the post-Newtonian 
and the tidal perturbation expansions used. This means 
that at the currently achieved expansion order, our ini- 
tial data are by no means guaranteed to be superior to 
the other approaches discussed above. Nevertheless, our 
approach can be systematically extended to higher order. 

Data that use analytic approximations valid on the 
entire hypersurface are difficult to find because such ap- 
proximations are usually valid only in certain regions of 
the hypersurface. One method to construct such data 
is by asymptotically matching two different physical ap- 
proximations in a region where they overlap. Asymp- 
totic matching, which for general relativity was developed 
in 0, |S HH ? consists of comparing the asymp- 

totic approximations of two adjacent solution inside of 
an overlap 4-volume. By comparing these asymptotic 
approximations, matching returns a map between coor- 
dinates and parameters local to different regions, which 
forces adjacent solutions to be asymptotic to each other 
inside the overlap region called the buffer zone. Since 
these solutions are close to each other in the buffer zone, 
it is possible to construct transition functions that merge 
these solutions, thus generating a smooth global metric. 

Alvi |3(j attempted to apply matching to binary sys- 
tems, but instead of matching he actually performed 
patching, because he set the physical approximations 



equal to each other at a 2-surface. As a result, Jansen 
and Briigmann [37^ found that Alvi's 4-metric, and in 
particular his extrinsic curvature, was hard to smooth 
with transition functions in the buffer zone, since the so- 
lutions were not sufficiently close to each other. Jansen 
and Briigmann thus concluded that these large discon- 
tinuities in Alvi's data renders them impractical for nu- 
merical implementation. In paper 1, binary systems were 
studied once more, but this time true asymptotic match- 
ing was implemented, thus successfully generating solu- 
tions that approach each other in the buffer zone. It 
was possible to smooth these solutions and thus to cre- 
ate useful initial data. The approach we used in paper 1 
was to match a PN metric in harmonic gauge [23, valid 
far from both black holes but less than a gravitational 
wavelength from the center of mass (near zone), to a 
perturbed Schwarzschild metric in isotropic coordinates, 
valid close to the holes (inner zone). 

In this paper, we study the method developed in pa- 
per 1 in more detail and we investigate and improve the 
resulting data. This improvement is mainly due to us- 
ing different coordinates in the PN approximation. In- 
stead of employing harmonic gauge, we now use a post- 
Newtonian metric in ADMTT gauge [S III E3| We 
will see that this leads to better matching. The reason 
for this improvement is that the 4-metric in ADMTT co- 
ordinates can be brought into a form that is very close 
to Schwarzschild in isotropic coordinates near each black 
hole. Therefore, we obtain a much smoother match be- 
tween this near zone metric and the inner zone one, which 
is also given in isotropic coordinates (see Fig. [2J. The 
results presented in Sec. IVII indeed show that match- 
ing with the ADMTT PN near zone metric works bet- 
ter than with the harmonic PN near zone metric. In 
addition, the ADMTT PN near zone metric is very sim- 
ilar to the perturbed Schwarzschild metric even near the 
black hole horizon, where PN approximations in principle 
break down. This similarity is due to the fact that, we 
use a resummed version of the ADMTT metric, obtained 
by adding higher order post-Newtonian terms. 

Another piece of evidence that suggests that matching 
works better between ADMTT and isotropic coordinates 
comes from looking at Hamiltonian and momentum con- 
straint violations. In Sec. IVllI we compare the Hamilto- 
nian and momentum constraint violations for the data 
of this paper and of paper 1. We find that the use of 
the ADMTT PN near zone metric leads to smaller con- 
straint violations than in paper 1. With the help of a 
transition function we can also obtain a global smooth 
4-metric which has errors on the order of the error intro- 
duced by the more accurate of the two approximations 
we match. 

By computing the constraint violations, we discover 
that both the data set of this paper and that presented 
in paper 1 might not be easily implemented in numerical 
simulations. In particular, we observe that the momen- 
tum constraint diverges near the apparent horizon of each 
black hole. This makes it difficult or impossible to ex- 
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cise a region inside each black hole, which contains all 
points where divergences occur. This divergence arises 
due to the choice of coordinates in the inner zone. Since 
we use the t = const slices of isotropic coordinates as 
spatial slices, the lapse goes through zero near the hori- 
zon. This zero in the lapse in turn leads to a blowup 
of some components of the extrinsic curvature and the 
momentum constraint. In order to circumvent this prob- 
lem, we construct a map from isotropic coordinates to 
new horizon-penetrating coordinates. A similar idea was 
presented in Ref. |4l|. but here we extend those ideas 
and provide explicit expressions for the transformed in- 
ner zone metric in a ready-to-implemcnt form. These 
new coordinates are identical Kerr-Schild coordinates in- 
side the horizon and become isotropic in the buffer zone. 
When this transformation is applied to either the data 
set presented in this paper or that presented in paper 1, 
the lapse of the new t = const slices is positive through 
the horizon, eliminating the divergence in the extrinsic 
curvature. This new coordinate system will make exci- 
sion easier, without changing the physical content of the 
data. 

The paper is organized as follows. Sec. ITTl explains how 
the spacetime is divided into zones and how matching will 
be implemented. Sec. IIIII describes the near zone post- 
Newtonian metric, while Sec. IIVI focuses on the inner 
zone metric. Sec. [3 performs the matching and provides 
a map between coordinates and parameters local to the 
near and inner zones. Sec. IVII gives explicit formulas 
for the global metric, introduces the smoothing functions 
and also decomposes this 4-metric into initial data for 
numerical relativity. Sec. IVIII compares the constraint 
violations of the initial data presented in this paper to 
those of paper 1. In Sec. I Villi we present an additional 
coordinate transformation which can be used to construct 
horizon penetrating coordinates. Sec. IIXI concludes and 
points to future work. Throughout we use geometrized 
units, where G = c = 1. 



II. DIVISION OF SPACETIME INTO ZONES 
AND MATCHING IN GR 

Consider a binary black hole spacetime and divide it 
into 4 zones. First, there are the so called inner zones 
Ci and C2 close to each black hole, where we can use 
black hole perturbation theory to obtain an approximate 
solution to the Einstein equations. These solutions are 
obtained under the assumption that the black holes are 
separated far enough that the influence of black hole 2 
is only a small perturbation near black hole 1. Second, 
there is the near zone C3 where PN theory should provide 
a good approximation as long as the black holes do not 
move too fast. Finally C4 denotes the far zone where re- 
tardation effects matter. These zones are shown in Fig.Q] 
and we summarize them in Table I. The quantities in- 
troduced in this table are defined as follows: fj n and r ou t 
are the approximate inner and outer boundary of each 





Far Zone (C 4 ) 


Near 

/ /^nner Zone BH 2 (C 2 ) 


Zone7c5~~~~~-~^ 

Buffer >v 
Zone(O l:j ) \ 

BH1 \ 






Inner Zone BH 1 (C\ / 

All 11 1. 1 I'll 1 1 V^l / / 


\ '-^A> J- * Buffer 

\ HH 1 Zone 







b 



FIG. 1: Schematic diagram of the near zone (dark gray), 
inner zones (light gray) and buffer zones (checkered). The 
near zone extends up to the inner boundary of the buffer 
zone, while the inner zone extends from the outer boundary 
of the buffer zone up to the black hole. Therefore, the buffer 
zones are spherical shells where the near zone and inner zones 
overlap. BH 1 and BH 2 (black) are separated by distance b. 
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Table I: Description of the division of the spacetime into 
zones |36il . 



zone respectively, as measured from the Ath black hole 
(A = 1 or 2) in the coordinate system local to it; e(„) 
is the expansion parameter used to find an approximate 
solution to the Einstein equations in that zone; ttla is the 
mass parameter of the A%\t black hole in the near and far 
zones; r a and f a are the radial distances as measured 
from the Ath. hole in the near and inner zones respec- 
tively; A is the typical gravitational wavelength; and b is 
the coordinate separation of the black holes. 

Inside of each zone we approximate the 4-metric with 
an monovariate expansion in £(«), which is defined in 
some coordinate system, and which depends on param- 
eters, like the mass of the black holes or the orbital 
angular velocity. The inner zones are both equipped 
with isotropic coordinates and have parameters (Ma, 
where the approximate solutions to the Einstein equa- 
tions are given by a perturbed Schwarzschild metric. On 
the other hand, the near zone is equipped with ADMTT 
coordinates (ADMTT gauge), where the approximate so- 
lution to the Einstein equations is a PN approximation 
with parameters uia and to. This PN approximation is 
valid in both C3 and C4 up to the order treated in this 
paper. 

As explained in paper 1, asymptotic matching consists 
of comparing the asymptotic expansions of adjacent ap- 



4 



proximate metrics in an overlap region or buffer zone. 
These buffer zones are defined by the intersection of the 
regions of validity of two adjacent metrics. The 2 buffer 
zones are delimited by the intersection of the inner zones 
C\ and C2 (of BHs 1 and 2), and the near zone C3. These 
buffer zones O12 = Ci H C2 and O13 = C\ (~l C3 can only be 
defined in the asymptotic sense [42j via wia « r/i « i 
and are the regions in which asymptotic matching will 
be performed. Formally, there is a third buffer zone 
(O34), given by the intersection of the near and far zone, 
but since the ADMTT PN metric at the order consid- 
ered here is valid in both C3 and C4, it is unnecessary 
to perform any matching there. The asymptotic expan- 
sions of the approximate metrics in the buffer zones Oaz 
are bivariate, because they will depend on 2 independent 
parameters: m^/r^ and r A /b. Thus, the errors in the 
asymptotic expansions valid in the buffer zones will be 
denoted as 0(p,q), which stands for errors of 0(m/b) p 
or errors of 0(rvt/&) 9 . 

Asymptotic matching then produces a coordinate and 
parameter mapping between adjacent regions. In this 
paper, these maps will consist of a transformation be- 
tween isotropic and ADMTT coordinates, as well as a 
set of conditions to relate Ma and f2 to tua and u>. With 
these maps, we can construct a piece- wise global met- 
ric by choosing a 3-surface inside the buffer zone, where 
we join the approximate solutions to the Einstein equa- 
tions together. Small discontinuities will be present in 
the global metric at the chosen 3-surface, due to its in- 
herent piecewise nature, but they are of the same order 
as the errors in the physical approximations, and thus 
controllable by the size of the perturbation parameters. 
Furthermore, if such parameters are chosen sufficiently 
small, these errors could, in particular, become smaller 
than numerical discretization errors. Ultimately, a C°° 
metric is sought after, so, since these discontinuities are 
small, the global metric can be smoothed by introducing 
transition functions. 



III. NEAR ZONE: AN ADMTT 
POST-NEWTONIAN METRIC 

In this section, we present the near zone metric in 
ADMTT gauge valid in C3 and C4 and expand it in the 
buffer zone O13. Analyzing matching in O13 will suffice 
due to the symmetry of the problem. The coordinate- 
parameter map for the overlap region O23 will later be 
given by a simple symmetry transformation. 

The ADMTT PN metric is obtained by solving the 
equations of motion for a binary black hole system with 
Hamiltonian dynamics |38j . The Hamiltonian is ex- 
panded in a slow-motion, weak-field approximation (PN) 
in ADMTT coordinates. This approach is similar to the 
standard Lagrangian PN expansion implemented in pa- 
per 1 and there exists a mapping between the two meth- 
ods [2^|. However, they differ in that the Hamiltonian 
formulation introduces a 3 + 1 decomposition of the 4- 



metric from the start, and thus it provides PN expres- 
sions for the lapse, shift, 3-metric and the conjugate mo- 
mentum. 

We follow the conventions of Ref. [3^ where the black 
hole trajectories from the center of mass of the system is 
given by 

= ^-b(n un = -^b(t"). (i) 

m m 

In Eq. l[TJl. m,A and m = TOi+m2 stand for the mass of BH 
A and the combined mass of the system in the ADMTT 
gauge, whereas the separation vector b(t) is given by 



b(t") = b {cos Lut", sin Lut",0) . 



(2) 



Note that Eq. (J2J assumes the black holes are initially in 
a circular orbit, which is a sensible approximation in the 
late stages of inspiral, since gravitational radiation will 
have circularized the orbit. The angular velocity of this 
orbit is given by Eq. (60) in Ref. |17|: 



w = y m/6 3 



1 



rn 



(3) 



with errors of 0(m/b) 5 ^ 2 . In Eq. M = miW2/m is 
the reduced mass of the system and b is the norm of the 
separation vector. This angular velocity is valid in the 
ADMTT gauge, which is different from the PN velocity 
used in paper 1, since the latter is valid only in harmonic 
coordinates. 

We further define the radial vector pointing from either 
black hole to a field point by 



(4) 



where the double primed variables (t" , x", y" , z") are in- 
ertial ADMTT coordinates measured from the center of 
mass. With these definitions, it is clear that f\" = f^'+b, 
where b points from black hole 1 to 2. We also introduce 
the unit vector 



n A " = r A "/r A ' 



(5) 



and the radial vector pointing from the center of mass to 
afield point r" = (x",y",z"). 

The particle's velocity vector is given by 

v A " = d t r A " = (-l) A+1 — bus (-smut", coswt",0). 
tcia 

(6) 

The post-Newtonian near zone metric in inertial 
ADMTT coordinates (Eqs. (5. 4)-(5. 6) in Ref. [H) can 
then be written as 
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A = ^ (3)i ", 



(7) 
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where we introduced a post-Newtonian conformal factor 



* = 1 



mi 
2?( 



ni2 



(8) 



and where the post-Newtonian lapse and shift (Eqs.(5.4)- 
(5.6) in Ref. |2g) are written as 



,(3) = 



(9) 



and 
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(10) 



These expressions are accurate up to errors of 0(m/rA) 2 - 
The 3-metric is conformally flat and in the vicinity of each 
black hole (near ta — 0) the metric is very similar to the 
metric of a black hole in isotropic coordinates. Note that 
Ref. uses units where 16irG — 1, but in Eq. Q we 
use units where G — 1 instead. 

Also note that the PN expressions above are not pure 
Taylor expansions in m/r or v/c. Instead, certain higher 
order PN terms have been added to make the metric sim- 
ilar to the Schwarzschild metric in isotropic coordinates 
near each black hole. This resummation formally does 
not change the PN accuracy of the metric, but in prac- 
tice it improves the PN metric since it yields a metric 
that is very similar to the inner zone black hole metric. 
The PN metric written in this form above even has ap- 
parent horizons near ta = ttla/^- Later when we plot 
our results we will use this resummed form of the metric. 
The reader may worry that there is a certain arbitrari- 
ness in adding higher order terms. However, we have used 
the following two criteria to minimize this arbitrariness: 
i) We only add terms that would actually be present at 
higher PN order (compare EHH |H E3 ) . h) The terms 
added have coefficients of order unity, and thus do not 
affect the usual PN order counting. Notice that point 
i) of course only means that we have added some higher 
order terms, not all higher order terms. 

The lapse and the spatial metric are similar to what 
was used in paper 1, when we expand them in the buffer 
zone. However, the shift given here is different from that 
in paper 1 and, thus, the maps returned by asymptotic 
matching will also be different. 

Now that we have the 4-metric in the near zone it is 
convenient to remove the frame rotation by performing 
the transformation 



t" = t, 

x = x cos cot — y sin cut, 

y = x smut + y cos Lot, 

z = z, 



where unprimed symbols stand for corotating ADMTT 
coordinates and double primed symbols correspond to 
the inertial ADMTT frame. After doing this coordinate 
change, the form of the near zone metric J7J) as well as 
the lapse and conformal factor remain unchanged, with 
only the shift picking up an additional lo x r term and 
now becomes 



m 1 
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(12) 



where va and ua are now time independent and equal to 
the double primed versions at t = 0. In these coordinates, 
the magnitude of the radial vectors become 
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rA 



2\ 1/2 



(13) 



where x\ = x — bm2/m and x% = x + bmi/m. 

Let us now make one final coordinate transformation 
where we will shift the origin to the center of BH 1, i.e., 



x = X\ = x 
.1 



V 



V, 



m 2 b 

m ' 

z = z, 



t' = t, 



(14) 



and r' A = ta- The single-primed coordinates then stand 
for the shifted corotating ADMTT coordinate system. 
This coordinate change is useful because we will later 
concentrate on matching in On. We will usually make 
explicit reference to x\ instead of x' to remind us that 
this distance should be measured from BH 1. This shift 
of the origin was not performed in paper 1 and, thus, 
the coordinate transformations that we will find will look 
slightly different. However, to 0(m/b) this transforma- 
tion should reduce to the results of paper 1 when the 
shift is undone. 

We now concentrate on the overlap region (buffer zone) 
013 . Inside O13 we can expand all terms proportional to 
l/r2 in Legendre polynomials of the form 
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(11) Substituting Eq. (|T5|) into Eq. J7J, we obtain 
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(16) 



where all errors are of order 0(2, 3) and where mi <C 
n <C 6. Eq. ((TBI) is denoted with a tilde because it is the 
asymptotic expansion in the buffer zone around BH 1 of 
the ADMTT metric. This metric has now been expressed 
entirely in terms of the shifted corotating ADMTT coor- 
dinate system (x 1 ', y',z',t'). 

Note that Eq. (|16D is a bivariate series, since it de- 
pends on two independent expansion parameters mi/r[ 
and r'l/b. Since both expansion parameters must be small 
independently, Ea. i|16|) is valid in a 4- volume defined by 
the size of the buffer zone. 



IV. INNER ZONE: A BLACK HOLE 
PERTURB ATIVE METRIC 

In this section we discuss the metric in the inner zone 
C\ of BH 1 and its asymptotic expansion in the overlap 



region O13. Since this metric will be the same as that 
used in the inner zone in paper 1, we will minimize its 
discussion and mostly refer to that paper. However, since 
we are using different notation than that used in paper 
1, we will summarize here the principal formulas. 



In the inner zone 1, let us use inertial isotropic coordi- 
nates, labeled by x 1 ' = (x',y',z',t'), which are centered 
at BH 1. Note that these coordinate are identical to the 
shifted corotated harmonic coordinates used in the near 
zone to zeroth order. Let us further denote the inner 
zone metric by g^pr, which will be given by a tidally 
perturbed Schwarzschild solution as 
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4ff 



(17) 



where fi = (x 12 + y' 2 + z' 2 ) 1 / 2 . This equation is identical 
to that used in the inner zone of paper 1 or Eq. (3.23) 
of Ref . [3(| . This metric was first computed by Alvi j3|j 
by linearly superposing a Schwarzschild metric to a tidal 
perturbation near BH 1. This perturbation is computed 
as an expansion in ei j2 = f' 12 /b and it represents the tidal 
effects of the external universe on BH 1. Since a simple 
linear superposition would not solve the Einstein equa- 
tions, multiplicative scalar functions are introduced into 
the metric perturbation, which in turn are determined 
by solving the linearized Einstein equations. With these 
scalar functions, the metric then solves the linearized 
Einstein equations and it represents a tidally perturbed 
Schwarzschild black hole. This metrics has been found 
to be isomorphic to that computed by Poisson |43| in 
advanced Eddington-Finkelstcin coordinates. 

Asymptotic matching will be easier when performed 



between metrics in similar coordinate systems. We, thus, 
choose to make a coordinate transformation to corotating 
isotropic coordinates x l = (x,y,z,i). The inner zone 
metric in corotating isotropic coordinates will be denoted 

by 9$ and is s iven b y 



■^00 

■ J oi 
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8jjH s i — H S 2 ^ 2 



(18) 



where ejjj, is the standard Levi-Civita symbol with con- 
vention ej23 = 1 and where 5^ is the Kronecker delta. In 
Eq. H18|) we use the shorthand 



H s t — 

H s \ = 

H S 2 = 

H t = 
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(19) 



r 



where f\ = f[ and the errors are still of 0(fi/6) 3 . These We now need to asymptotically expand the inner zone 
shorthands for the different components of the metric are metric in the buffer zone, i.e., as Mi/fi <C 1. Doing so 
identical to those used in paper 1. 
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we obtain 



g(l) 
»01 



~(D 
y 02 



-(1) 
^03 







2M, 2m 2 _ 2 

-1 + — + -73-fl-P2 , 

ri 6 J \ri 
+fi 2 (x 2 + y 2 ), 

2m 2 /m , 2 _^ / 2Mi 



2m 2 /m 
6~ 

2Mi , 2m 2 _ 2 



63 



1 



ri 



63 ^ 2 1 - 



The asymptotic expansion of gj^ will be denoted by 



(20) 



g^p. Note that this asymptotic expansion is a bivariate 
expansion in both f±/b and M\jf\. In other words, it 
is the asymptotic expansion in the buffer zone to the 
approximate solution in the inner zone. 



V. MATCHING CONDITIONS AND 
COORDINATE TRANSFORMATIONS 



In this section, we find the coordinate and parameter 
maps that relate adjacent solutions. Since the coordinate 
systems are similar to each other in the buffer zone, we 
make the following ansatz for the transformation between 
coordinates and parameters 



t 

Mi 
CI 



i + (^) i/2 71 (^) + (^) 72 (^) + (^) 3/2 73 (^) 

i + (^) 1/2 Cl (^) + (^) C2 (^) + (^) 3/2 C3(^) 



x 3/2 

) T 3 (^ ) 



/m 2 \ 1 /2 m 2 (m 2 \ 



3/2 



6 

v 6 ; 



/m 2 \!/2 TO2 /m 2 \3/2 

I f ( — ! Kl + — K 2 + y — 



K3 



(21) 



where, as in paper 1, Xi,2,3, 71,2,3, Ci,2,3 and 71,2,3 are 
functions of the coordinates, whereas 771,2,3 and Ki, 2j3 are 
coordinate independent. Note that this ansatz is slightly 
different from the one made in paper 1, because here we 
have shifted the origin of the near zone coordinate sys- 
tem, so that to zeroth order it agrees with the coordinates 
used in the inner zone. In order to determine these maps, 
we enforce the matching condition, i.e. 



dxP dx 17 
dx^ 1 dx v ' ' 



(22) 



which leads to a system of first-order coupled partial dif- 
ferential equations, which we solve order by order. 

The system to 0(m/b) provides no information be- 
cause both metrics are asymptotic to Minkowski space- 
time to lowest order. The first non-trivial matching, oc- 
curs at 0(m/6) 1 / 2 . The differential system at this order 
and at order 0{m/b) are similar to those obtain in pa- 
per 1 and we, thus, omit them here. The solution up to 



0(m/b) is given by 



to 2 / x'\ /m 2 \3/2 , ' 
1 + X ^26 +(-) X3(X ] 



1 + (it) \J 



m 2 m 2 , 
+ a/ ~6~V^ ' 



™2\3/2 / 

-y) 73 {x" ) 



X 



(23) 



with 771 = and m = 0. We should note that in solv- 
ing the differential system to 0(to/6) 1 / 2 and 0(m/b) we 
have explicitly required that the coordinate system be 
not boosted or rotating with respect to each other. Fur- 
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thermore, we found the particular solution where the 
t' = and the i — slices coincide, since t 1 = im- 
plies i — 0. We should also note that Eq. Ij23(l of this 
paper is different from the transformation found in pa- 
per 1 because in this paper we have shifted the origin 
of the near zone coordinate system, which simplifies the 
transformation. 

Matching at 0{m/b) 3 / 2 needs to be performed only 
on the Oi components of the metric to get the first non- 



trivial correction to the extrinsic curvature. This order 
counting does not follow the standard post-Newtonian 
scheme and it is explained in more detail in paper 1. 
Since the near zone shift is different from that used in 
paper 1 , the differential system obtained by imposing the 
matching condition of Eq. I|22[) is also different. This 
system is given by 



t'T 3 . x , - x'xs,t' 
t'T 3 , y ' - y'lz,t> 



t'T 3<z < - z'C 



3,t' 



2/V 

b 2 



m 
m 2 
rn 
m 2 
7m i 
2m 
m y'z' 
m 2 b 2 



"mi / 


'b 3 


2m ' 


y 3 


7m i 


3 


4m 


~ 2 


- Ua 


- 1 



m 2 
rn 

m 2 
m 



b 2 
7m 1 
2m 





y' 


" m 2 






m 



mi 
2m 



mi 
2m 
3 

+ 2 



'l 

m 2 
m 



K 2 "+ 

m 2 
rn 



1 - 

r 

~ 2 



mi 
2 m 



7 'mi 
2m 



b 2 



x 



JL 

2m 



m 2 



(24) 



Note that k 2 represents the degree of freedom associated 
with the angular velocity. For simplicity, we now set n 2 = 
0. Also note that rj 2l i] 3 and k 3 never enter the matching 
conditions to this order, so we set those quantities to 
zero. All the matching relations between the parameters 
local to different approximation have now been specified 
to 0(m/b) 2 , i.e. 

flttuj, Mi w mi. (25) 



With these matching relations, we can now find a par- 
ticular solution to Eqs. (|24|l and this set will produce 
matching. Once more, we choose constants of integra- 
tion such that the t' = and the t = slices coincide. 
Doing so we obtain the solution 




where errors are of (9(2,3). replaced r[ by 

As in paper 1, this coordinate transformation is singu- , 

lar at r[ = 0. Note, however, that in Eq. (gEJ we have £i = y r i + 6™ 2 - (27) 
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This replacement amounts to adding higher order terms 
to the coordinate transformation, which have no effect in 
the buffer zone at the current level of accuracy. Yet it 
has the advantage that the resulting coordinate transfor- 
mation is now regular at r[ — 0. Thus, the replacement 
does not affect asymptotic matching (in the buffer zone) . 
It merely leads to an inner zone metric expressed in a 
better coordinate system. 

In paper 1 there was one term in the coordinate trans- 
formation that was singular at the location of the holes, 
whereas here there are 3 such singular terms. This in- 
crease in singular behavior is due to the ADMTT shift 
having more poles than the shift of paper 1. The shift 
of the origin into the complex plane, as introduced in 
Eq. 1)27(1. removes this singular behavior from the real 
line, while only introducing uncontrolled remainders at a 
higher order. We should also note that in Eq. I|26|l . we 
undid the shift given by Eq. I|14|) . so that we measure 
distances from the center of mass of the system. 

The coordinate transformations presented in Eq. 12(j|) 
is only valid in O13, but we can find the transforma- 
tion in 023 by a simple symmetry transformation. This 
transformation is given by the following rules: substitute 
1 <-> 2 and 



y 



-y, 



(28) 



The set of matching conditions plus the coordinate 
transformation presented in Eqs. (|25|l and i|26|) are the 
end result of the matching scheme. Now that we have 
an approximate global metric, we can decompose it to 
provide initial data for numerical simulations. This de- 
composition consists of a foliation of the manifold with 
spacelike hypersurfaces, where we choose the slicing given 
by t = = i . Notice that we have found a coordinate 
transformation that is consistent both with this foliation 
and with the condition rriA <C rA *C valid in the buffer 
zone. 

Note that there is some freedom in the matching 
scheme at the order of accuracy of this paper. This free- 
dom is rooted in that a different choice for the matching 
parameters (M and f2) could be made. However, except 
for the K2 term in f2, this would have no effect on the 
coordinate transformation or the 4-metric at the orders 
considered here. Furthermore, one can show that the 
value of K2 does not affect the physics of the initial data 
we will construct below. Above, the choice of K2 = was 
made in order to simplify the relation between the pa- 
rameters of the two approximations. Specifically, it was 
chosen such that the angular velocity parameters of the 
two approximations are equal. This is, however, not the 
only possible choice. Fortunately, the choice of K2 will not 
affect the initial data we will construct later, since both 
the 3-metric and the extrinsic curvature will change only 
by uncontrolled remainders of higher order. This hap- 
pens for the following reasons. A non-zero K2 would only 
change the spatial part of the coordinate transformation 
at 0(m/6) 3 / 2 , and thus would have no effect on the slic- 
ing. Since the leading order non-zero piece in the extrin- 



sic curvature is already of 0(m/&) 3 / 2 and proportional to 
fi, any change due to spatial coordinate transformations 
or parameter changes at 0(m/b) 3 ^ 2 will cause changes 
only at even higher order. For the 3-metric, coordinate 
transformations or parameter changes at 0(m/b) 3 ^ 2 will 
make a difference at order 0(m/b) 3 ^ 2 in the metric. We 
can, however, again neglect these changes, since the 3- 
metric was only matched up to 0(m/b). In conclusion, 
this means that any allowed change in the matching pa- 
rameters will not affect the physics of the initial data sets 
constructed below. 



VI. CONSTRUCTING A GLOBAL METRIC 

In this section, we present the global metric, by per- 
forming the coordinate transformation found in the pre- 
vious section on the inner zone metric. The piecewise 
metric in corotating ADMTT coordinates is given by 



(global) 
1 \iu 



o (1) 


in 




(2) 


in 


c 2 


(3) 

9f.v j 


in 


c 3 



(29) 



In Eq. (|29|l . gj?J denotes the near zone metric given in 
Eq. Q, but in unprimed corotating ADMTT coordi- 
nates, whereas gj^J and gffi stand for the inner zone met- 
rics of BH 1 and 2 transformed to corotating ADMTT 
coordinates via Eq. (|26|l . Explicitly, the inner zone met- 
ric is given by 



„(1,2) 



(1.2) jfi jv 



where Jj? is the Jacobian of Eq. J5BJl, namely 



JS 



dx^ 
ft? 7 ' 



(30) 



(31) 



In Eq. (|3Ufl . g~ is the inner zone metric in corotating 



isotropic coordinates given by Eqs. (|18fl. while g^J can 
be obtained by applying Eq. QgJ to gj£. Eq. (J20J is 
also expanded in paper 1 with the substitution fify — 
g9p. Note that the Jacobian given in Appendix lAl will 
be different from that presented in paper 1 because the 
coordinate transformation is different. 

In the following we present plots of different metric 
components at time i = t' = 0. These plots are represen- 
tative samples of the initial data and will show how well 
the matching procedure works. In most of these figures 
we choose to plot along the axis that connects the holes 
(the x-axis). We have checked that the behavior along 
other axis is qualitatively similar to that along the x-axis, 
as evidenced by the contour plots included in the section. 
Furthermore, we plot only the relevant components be- 
cause either the other components vanish along the x-axis 
or they present similar behavior. Thus the components 
and axis chosen present the general behavior of the data. 
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For these plots, we choose two different physical sys- 
tems. The first system consists of equal mass black 
holes separated by b = 10m. The two black holes are 
located at x = ±5m and are surrounded by apparent 
horizons, which approximately are coordinate spheres of 
radius m/4. For this system, some quantities, like the xx- 
component of the metric, are reflection symmetric under 
x — > —x and we then plot only the positive x-axis, omit- 
ting the inner zone 2 data. The second system consists 
of two holes with mass ratio m^/mx = 0.25 separated 
by b = 20m. Here the holes are located at x = 4m 
and x = —16m, and surrounded by apparent horizons 
with radii of approximately mi/2 and 7712/2 for holes 
1 and 2 respectively. In general, a dashed line corre- 
sponds to the post-Newtonian near zone metric, while the 
dotted or dot-dot-dashed lines represent the perturbed 
Schwarzschild inner zone metrics. 

Most figures will contain error bars, which in the case 
of the metric are given by 







5mim 2 






U " 


-f) 


4 b 


in' 





mim 2 



S 2 ' 

Einner,A = 2— ; (32) 

where S = r*i + r 2 + b. These error bars were obtained by 
analyzing the next order term in the ADMTT PN met- 
ric and in the black hole perturbation approximation. In 
particular, the PN error is different from that used in 
paper 1. Here the first term comes from the expansion of 
the conformal factor and the last two from the 2 PN part 
of the trans verse-traceless term in the ADMTT metric. 
The full PN error, given by the 3.5 ADMTT PN metric, 
is much more complex than the simplified approximation 
used in E near , where the former contains some terms 
that scale for example as m 2 A /r\ and m\/{Sr A ). How- 
ever, here we have neglected those terms and we have 
checked that E near models well both the magnitude and 
functional behavior of the full error in the regions plot- 
ted in all figures with fractional errors of less than 1% 
everywhere. 

In Fig. |21 we plot the xx-component of the piece- wise 
metric along the positive x-axis. Observe that the in- 
ner zone metric of black hole 1 (dotted line) is close to 
the post-Newtonian near zone metric (dashed line) in the 
buffer zone of black hole 1 given by mi <C f\ -C b. Re- 
call that the definition of the buffer zone in the theory of 
asymptotic matching is inherently imprecise in the sense 
that as one approaches either the inner or outer radius 
of the buffer zone shell, one of the two approximations 
(i.e. near zone or inner zone metric) has much larger er- 
rors than the other. In our case we need both fji/b and 
mA/rA to be small, which will not be the case at either 
end of the interval vha *C Ta ^ b. Hence we expect the 
best agreement between the two approximations to oc- 
cur somewhere in the middle of each buffer zone shell. In 
fact the best agreement should occur where each of the 
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FIG. 2: This figure shows the xx-component of the near 
zone (dashed line) and inner zone 1 metric (dotted line) along 
the positive x-axis in ADMTT coordinates. The black holes 
are located at x/m w ±5 and the perturbative parameter is 
m/b — 1/10. Observe that the metrics become asymptotic 
to each other in the buffer zone, whose approximate center is 
roughly 4m away from each black hole. Also plotted are esti- 
mates of the error in the approximations. Note that the off- 
diagonal components of the 3-metric vanish along the x-axis 
for symmetry reasons, and that the yy- and zz-components 
are very similar to the xx-component. 

two approximations have roughly the same error. From 
Fig. we see that the inner zone metric of each black 
hole agrees with the near zone metric within error bars 
in the middle of each of the two buffer zones, where both 
have comparable errors. This agreement means that we 
have performed successful matching. However, our re- 
sult in Fig. |3 is even better than this, since the near 
zone 3-metric is close to each of the inner zone 3-metrics 
even near the black holes, where post-Newtonian the- 
ory has large errors and is expected to fail. The reason 
for this somewhat surprising success near the black holes 
is that we have used resummed post-Newtonian expres- 
sions [Eqs. 0, ©, © and JTOJl] in the ADMTT gauge, 
which are very close to the isotropic Schwarzschild met- 
ric used in the inner zone. Note that the off-diagonal 
components of the 3-metric vanish along the x-axis for 
symmetry reasons, and that the yy- and zz-components 
are very similar to the xx component. 

The quality of the asymptotic matching (in the buffer 
zone) is equally good for all other components of the 4- 
metric. This is evidenced in Figs.[3Iand|Hwhich show the 
near and inner zone lapse and y-component of the shift, 
computed from the 4-metric (|29|l via 

P = g ij g oj . 

a = = (g ok p k - g 00 ) 1/2 . (33) 

As one can see from Fig. [3] the near (dashed line) and 
inner zone 1 lapse (dotted line) agree well in the buffer 
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FIG. 3: Plot of the lapse along the positive :r-axis for the 
inner zone and near zone metrics. Note that the lapses are 
similar to each other not only in the buffer zone, but also near 
the hole. 
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FIG. 4: Plot of the y-component of the shift along the positive 
a;-axis for the inner zone and near zone metrics. Observe that 
although the inner zone and near zone curves are similar in 
the buffer zones, in this case they differ near the black holes. 
Also note that the x- and ^-components of the shift vanish 
along the i-axis. 



FIG. 5: This figure shows the xx component of the piece- 
wise metric for a system of unequal masses and at larger sep- 
arations. In particular, we have used mi/m\ = 0.25 and 
b = 20m. Note that the matching improves since the pertur- 
bation parameter m/b is smaller. 



addition, for most components (except the shift) the re- 
summed post-Newtonian expressions [Eqs. J7J|, ©, © 
and i|l(Jfl ] we use in the near zone, follow the inner zone 
behavior even near the black holes, where this is not nec- 
essarily expected. 

We have also tested the coordinate and parameter 
maps for system with unequal masses and found that 
matching is still valid. Figure [S] shows the aix-component 
of the piece-wise metric along the entire £-axis for the 
case of m2/ 'mi = 0.25 and b — 20m. Observe that the 
matching improves in this case because the perturbation 
parameter m/b has become smaller. Also note that the 
general features of the matching are the same as those 
presented in the equal mass case. We have plotted along 
the entire x-axis and show both the metrics of inner zone 
1 (dotted line) and 2 (dot-dot-dashed line). 

Since we are interested in initial data on the i = t' = 
slice we also have to discuss the extrinsic curvature given 
by 



zone, located about Am away from each black hole. Fur- 
thermore, the resummed expression [Eq. ©] we use for 
the near zone lapse also is very close to the inner zone 
lapse near each black hole, where post-Newtonian theory 
has large error bars. In Fig. 0] we see that the near zone 
shift (dashed line) is close to the inner zone 1 shift (dotted 
line) in the buffer zone as expected after matching there. 
However, in the inner zone the post-Newtonian near zone 
shift is not valid anymore and deviates strongly from the 
inner zone shift. 

Thus in summary asymptotic matching works well for 
all components of the 4-metric in the buffer zone. In 



Kij = (d t g i:i ~ £ gij) . (34) 

In the inner zones we use this equation to numerically 
compute Kj-j and '. Note however, that the inner 
zone extrinsic curvature in corotating isotropic coordi- 

(1) ICC 

nates can also be found in paper 1, where = 

KSV in the new notation. 

v 

The near zone extrinsic curvature can easily be com- 
puted analytically from the near zone post-Newtonian 
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FIG. 6: This figure shows the xy component of the near 
zone (dashed line) extrinsic curvature, as well as the inner 
zone curvature (dotted line) obtained via black hole pertur- 
bation theory. The plot is for equal mass black holes with 
m/b = 1/10. In the buffer zones, where the near and inner 
zone approximations have comparable error bars, the two ap- 
proximations agree within error bars. All other components 
vanish along the a;-axis. 



FIG. 7: This figure shows the xy component of the piece-wise 
extrinsic curvature for a system of unequal masses (mz/mi = 
0.25) at separation b = 20m. Note that the scale of the y-axis 
has been reduced. The curves are closer to each other in this 
case because the perturbation parameter m/b has decreased 
by the choice of parameters. 



4-metric in ADMTT coordinates. The result [l7| 



K 



(3) 



* 2 S 2r 2 - 



A=l 



PA n A + PA n A 



PA n A^mn (5 lj - n A n A ) 



(35) 



is of Bowen-York form, and v A and n l A denote the particle 
velocities and directional vectors already introduced in 
Sec. IIIII On the t' = slice these quantities become 



2 m 2i 

V 1 = OJ 0, 

m 



-u^b, v A = v\ = 0, (36) 
m 



and 



e A 



x 



TA 

i = o. 



m 



mi 

b, 

m 



(37) 



In Fig.[S]we plot the xy component of the extrinsic cur- 
vature along the positive x-axis in the equal mass case. 
For reasons of symmetry all other components vanish 
along the cc-axis. Observe that the near zone solution 
(dashed line) diverges faster than the inner zone solution 
(dotted line). From this figure we see that even though 
the resummed ADMTT post-Newtonian expansion mod- 
els well the 3-metric of a perturbed Schwarzschild hole, 
the ADMTT and inner extrinsic curvatures do not agree 
well in the inner zone, just like in the case of the shift. 
However, in the buffer zone where both inner and near 
zone results have similar errors, both approximations 
roughly agree. On the side away from the companion 



black hole this agreement is within error bars, while be- 
tween the holes the two error bars do not quite overlap. 
This suggests that a separation of b = 10m is at the 
edge of validity of the matching scheme employed here. 
Notice, however, that the error bars 



Erpn 



Erbhpt,. 



TUT- 



5 mim 2 ( 1 
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(38) 



plotted here (for v — (m/b) 1 / 2 ) are only approximations 
of the true errors. These approximate errors were ob- 
tained by differentiating the error bars for the metric in 
Eq. (|32|l with respect to time. Note that the absolute er- 
ror in the extrinsic curvature is 0(m/b)^' 2 smaller than 
that of the metric because the first non-trivial term in 
the extrinsic curvature appears at 0(m/b) 3 ^ 2 . 

Similar results for the extrinsic curvature are also ob- 
tained for different separations and mass ratios. Figure[7| 
shows the xy component of the piece-wise extrinsic cur- 
vature along the x-axis for the case of ma/rai = 0.25 and 
b = 20m. The near zone curvature (dashed line) agrees 
within error bars with the inner zone 1 curvature (dot- 
ted line) and the inner zone 2 curvature (dot-dot-dashed 
line) in buffer zones 1 and 2 respectively. Observe how- 
ever, that the agreement worsens outside the buffer zone, 
e.g. as we approach either hole. 
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Transition Functions 

In this subsection, we construct transition functions 
that allow us to remove the piece- wise nature of Eq. i|29|) 
and merge the approximate solutions smoothly in the 
buffer zone. This smoothing is performed by taking 
weighted averages of the two approximations. This pro- 
cedure is justified by the fact that both approximations 
are equal to each other in the buffer zone up to uncon- 
trolled remainders (asymptotic matching theorem 
which were neglected in the approximations used. 

In the following we assume that the middle of each 
buffer zone is located around 



m^rru^O.Sm, b=10m 



< = (b 4 m 2 A /m) 



1/5 



(39) 



which corresponds to the distance from the black holes 
where the error bars of the adjacent approximations are 
comparable. 

The transition functions we use are all of the form 



f(r) = 



0, r < r 
i {l + tanh[f (tan(^(r 



1, 



tan(^(r— r )) 

r > ro + w 



)]}■ 



ro)) 

r < r < r + w 



(40) 



which is a function which smoothly transitions from zero 
to one in the region ro < r < ro+w. The parameters used 
here are as follows: ro defines where the transition begins; 
w gives the width of the transition window; q determines 
the point where the transition function is equal to 1/2 
(this happens at r 1 / 2 = ro + (2w/ir) arctang); the slope 
of the transition function at this point is s(l + q 2 )/(2w) 
and thus can by influenced by choosing s (and of course 
w). 

The global merged 4-metric is then given by 
gjgoW = G(x){F l (f 1 )g^ + [l-F 1 (f 1 )]g^} 
+ [l-G(x)]{F 2 (f 2 )g$ + [l 

-F2(r*)]9$}> ( 41 ) 
where we have introduced 

F A (f A ) = f(f A ), (42) 

with 

r = QArf, w = 3.5r^, q = 0.2, s = b/m, 

(43) 

and also 

G(x) = f{x), (44) 

but with 

b(rri2 — mi) b 
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2m 2 
q = 1, s = 2.5. 
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FIG. 8: This figure shows the transition functions used to 
merge the different pieces of the 4-metric. The transition 
function F\(f\) with the transition parameters of Eq. I43H 
is shown as a solid line, while the function G(x) with the 
parameters of Eq. 1451 is shown with a dashed line. The 
transition function Qi(fi) (dotted line) defined in Eqs. 1571 
and 1581 1 is used in Sec lVIlH to construct horizon penetrating 
coordinates via the coordinate transformation 15511 . 



The transition functions Pi(fi) and G{x) are shown in 
Fig. [SI for m A — to/2, b = IOto and the parameters of 
Eqs. |@3J and gSJ, respectively. 

The global lapse, shift and extrinsic curvature are 
merged using the same transition functions as for the 4- 
metric above. With these transition functions, the global 
metric is a mixture of inner zone and near zone metric 
in a transition region given by 0.4r^ J < f A < 3.9r^ f . On 
the other hand, the global metric becomes identical to 
the inner zone metric in the region f A < 0.4r^, while 
it becomes equal to the near zone metric in the region 
f A > 3.9r^. Note that although the transition function 
are identical to those used in paper 1, the transition pa- 
rameters chosen here are slightly different. This change 
is because the inner zone metric is very similar to the 
near zone metric close to the black holes in ADMTT co- 
ordinates, so that the transition region has been moved 
closer to the black holes. 

The exact choice of a transition function is to a certain 
degree arbitrary and could in principle be changed. How- 
ever, the resulting global metrics generated by any other 
reasonable transition function should look very similar, 
and in fact be identical up to uncontrolled higher order 
post-Newtonian and tidal perturbation terms. Also note 
that the functions presented in this section are general 
enough to perform a smooth merge for systems with dif- 
ferent masses and separations. 

In Fig. (J5J), we present the zx-component of the metric 
around BH 1. In this figure, the dotted curve corresponds 
to the inner zone metric, the dashed line is the near zone 
metric and the solid line is the merged global metric. The 
error bar in the global metric is given by the smallest of 
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FIG. 9: In this figure we show how the transition function 
takes the xx component of the near zone metric to the inner 
zone solution around BH 1. Observe that the transition func- 
tion is smooth and does not introduce kinks into the global 
solution. 



FIG. 11: In this figure we plot a system of unequal masses 
(milm\ = 0.25) and a larger separation (6 = 20m). We use 
the same parameters in the transition function as in the equal 
mass case. 
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FIG. 10: In this figure we show the x-component of the met- 
ric, divided by W*. The dashed line is the near zone metric, 
the dotted line is the inner zone metric and the black solid 
line is the global metric with transition functions. 



the error bars of the respective approximations. The pur- 
pose of this figure is to show explicitly that the transition 
function effectively merges the different approximations 
in the buffer zone, where the errors are comparable. 

Even though the near zone 3-metric models the inner 
zone 3-metric well near the black hole horizons, the two 
metrics diverge differently at ta — 0. In Fig. llUI we plot 
the xx-component of the metric divided by \I/ 4 . This de- 
nominator removes the divergence of the near zone met- 
ric, which now becomes identically unity, while showing 
that the inner zone metric differs by an amount equal to 
the size of the tidal perturbation near the horizon. The 
reason that the inner zone metric still diverges after di- 



viding by 'J 4 is due to the fact that the tidal perturbation 
piece of the inner zone metric is divergent at r a = 0. This 
happens because the tidal perturbation piece is valid only 
for finite space-like separations from the horizon, while 
ta = which is located at the inner asymptotically flat 
end inside the black hole, is infinitely far from the hori- 
zon. This means that the tidal perturbation piece of 
the inner zone metric is not valid near ta = 0. Yet, if 
this spurious tidal perturbation is removed, the result- 
ing 3-metric of the Schwarzschild background will again 
agree well with the resummed post-Newtonian near zone 
3-metric. 

The same analysis can be performed for systems of 
unequal masses and at different separations. In Fig. 1111 
we plot the xx component of the metric along the x- 
axis for a binary with mass ratio m-ij mi = 0.25 and 
at a separation b = 20m. The transition is now easier 
than before (note the scale of the y-axis), because the 
perturbation parameter has become smaller. Also note 
that we have used the same parameters in the transition 
function as those used in the equal mass case and the 
transition is still smooth. 

With these transition functions we can construct a 
smooth global metric, as shown in Fig. There are 
no discontinuous features in the global metric. Just as 
in the case of the metric, to generate a global extrinsic 
curvature we will use the same transition function of the 
previous section with the same parameters. Figure 1131 
shows the global extrinsic curvature with the transition 
functions. Note that we could have computed the extrin- 
sic curvature directly from Eq. (|41|t . but this curvature 
would have been very similar to the curvature merged 
with transition functions. The difference between these 
two curvatures lies in the derivatives of the transition 
functions, but these terms will be of the same order or 
smaller than the uncontrolled remainders in either ap- 
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isolines in xy-plane 




FIG. 12: Contour plot of the xx-component of the global 
metric in the xy-plane around black hole 1 for the equal mass 
system. Apart from the singularity at x = ±5m, y = the 
3-metric is everywhere smooth. The different line styles cor- 
respond to different isolines of constant metric value. 



K isolines in xy-plane 




FIG. 13: Contour plot of the xy component of the global 
extrinsic curvature in the xy-plane around hole 1 for the equal 
mass system. The different line styles correspond to different 
isolines of constant extrinsic curvature. 



proximation in the buffer zone because of the parameter 
choice in Eqs. and |@SJ. 

The 3-metric g^ lobal > an d the extrinsic curvature given 
in this section can now be used as initial data for 
black hole binaries. Recall, however, that these data 
are solutions to the Einstein equations only approxi- 
mately and, thus, they only approximately satisfy the 
constraints. An analytic estimate of the constraint vi- 
olation and physical accuracy of this data is given in 
Table II. Note that in the near zone, both the errors 
in the physics and the constraint violations are given 
by the neglected terms in the post-Newtonian approx- 
imation, which scale as 0(m/rA) 2 - On the other hand, 
in the inner zone, the leading error in the physics is 
due to the terms neglected in the perturbation, which 



Zone 


Constraint Viol. 


Physical Error 


Inner zone BH1 (Ci) 
Inner zone BH2 (C 2 ) 
Near zone (Cz) 


0[(m/&) 5/2 (fx/fa) 2 ] 
0[{m/bf' 2 {f. 2 /bf] 
0(m/r A ) 2 


0[(m/b)(n/b) 3 ] 
0[{m/b){f 2 /b) 3 ] 
0(m/r A ) 2 



Table II: This table shows an order of magnitude estimate of 
the constraint violation of this data, together with the error 
in the physics. 



scales as 0[(m/b)(r /b) 3 ]. The inner zone constraint vi- 
olations, however, are smaller, because the perturbed 
Schwarzschild metric used here, satisfies the Einstein 
Equations up to order 0[(m/b) 5 ^ 2 (r /b) 2 ]. 

In order to obtain initial data that satisfy the con- 
straints exactly, it will be necessary to project the data 
given in this paper to the constraint hypersurface. How- 
ever, since these data are already significantly close to 
this hypersurface, sensible projection methods should not 
alter much the astrophysical content of the initial data. 
Furthermore, as stressed earlier, if this constraint viola- 
tion is smaller than discretization error, these data could 
be evolved without any projection. 

There are a couple of caveats that need to be discussed 
in more detail. First, as in paper 1, asymptotic match- 
ing will break down when the separation of the bodies 
is small enough that the near zone disappears between 
the two black holes. An approximate criterion for when 
this happens, corresponds to the separation where the 
middles of the two buffer zone shells touch for the first 
time. If we define the middle of each buffer zone as in 
Eq. (|39|l . this happens when the two spheres of radius 
rf 1 and r?f centered on each of the black holes, touch 
for the first time (roughly b « 8m for an equal mass sys- 
tem). Furthermore, the tidal perturbation used in the 
inner zone metric is only valid for small spacelike sepa- 
ration from the horizon. This criterion implies that the 
tidal perturbation is good roughly for 



tua m rriA 
— 2b <<fA<< ^ b - 



(46) 



Therefore, numerical simulations will need to excise the 
data somewhere inside the apparent horizons of the black 
holes before evolving it, since the inner zone metric is 
not valid all the way up to ta = 0. Notice, however, 
that these limitations are due to the approximations and 
coordinates used and not due to asymptotic matching. 



VII. CONSTRAINT VIOLATIONS 

The initial data constructed here are based on approx- 
imate solutions and, thus, they do not solve the Einstein 
equations exactly. We have estimated that the largest 
error in the constraints of the full theory occurs in the 
buffer zone and is at most 0(iriA/b) 2 . This error can be 
sufficiently small compared to other sources of numerical 
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FIG. 14: This figure shows the violation of the Hamiltonian 
constraint for the equal mass system along the positive :r-axis 
for both the data presented in this paper (solid black line) 
and that of paper 1 (dot-dashed black line). 



FIG. 15: This figure shows the violation of the Hamiltonian 
constraint for the equal mass system along the positive a;-axis 
for the near zone (dashed line), the inner zone (dotted line) 
and the global (solid line) data with transition functions. 



error such that solving the constraints more accurately 
may not be required. On the other hand, the initial data 
could be used as input in one of the many conformal 
decompositions 14] to compute a numerical solution to 
the full constraints. In this section, we study the con- 
straint violations of the data presented in this paper and 
we compare them to what was obtained in paper 1. 

In Fig. El we show the Hamiltonian constraint vio- 
lation along the positive x-axis for both data sets. The 
solid line is the constraint violation produced by the data 
presented in this paper, while the dot-dashed line is that 
produced by the data of paper 1. Observe that far away 
from both holes, both data sets (paper 1 and this pa- 
per's) have a constraint violation that approaches zero. 
This feature arises because both global metrics are given 
by a 1 PN solution in the near zone, which has a con- 
straint violation that decays when the expansion param- 
eter ttiaIta becomes smaller. Note, however, that the 
constraint violation in ADMTT gauge is much smaller 
than in harmonic gauge. The reason for this difference is 
that our resummed PN expressions satisfy the Hamilto- 
nian constraint up to errors of O^a/ta) 3 in ADMTT 
gauge, while the harmonic gauge expressions have errors 
of 0(m A /r A ) 2 . 

Also in Fig. El observe that close to the holes, the vi- 
olation is given by that of the inner zone solution, which 
is small near the horizon, but increases as we move away 
from the holes. The largest constraint violation occurs in 
the buffer zone, where we transition from inner to near 
zone expressions. Note that this maximum in the vio- 
lation is mainly caused by the transition functions. As 
can be seen in Fig. 1151 both the inner and near zone so- 
lutions individually have smaller constraint violations in 
the buffer zone than the global curve. This, however, 
is not an indication of a bad choice of transition func- 
tions. Rather the inner and near zone solutions have 
smaller constraint violations than expected, because we 



have kept some extra higher order terms. Recall that 
matching in the buffer zone was only performed up er- 
rors of 0{mA/rA) 2 ■ Hence the two solutions differ by 
this amount, and thus the error in the physics is of the 
same order as well. Nevertheless, the inner zone solution 
has constraint violations of C^m^/r^) 5 / 2 only, and the 
errors in the Hamiltonian constraint in the ADMTT near 
zone solution are only of 0(iriA/rA) 3 ■ When these two 
solutions (which are equal to up to order 0(m,A/rA)) are 
averaged with a transition function we obtain a new solu- 
tion which differs from both inner and near zone solution 
by O^tia/ta) 2 ■ Therefore, we expect the error in the 
constraints of the averaged solution to be of O^a/^a) 2 
as well, and thus larger than the error in the individual 
solutions. 

From Fig. El we see that the Hamiltonian constraint 
for the data of this paper is smaller than the data of 
paper 1. This decrease in the constraint violation is an 
indication that the matching performed in this paper pro- 
duces near and inner zone solutions that are closer to 
each other in the buffer zone. Thus, the transition func- 
tion has to do less work to join the solutions together, 
therefore introducing less of a constraint violation. We 
should point out that while the functional form of the 
transition function in paper 1 and in this paper is identi- 
cal, the parameters used are different. The biggest differ- 
ence is that here we use a smaller transition window w. 
Thus stronger artificial gradients and larger constraint 
violations might be expected in this work. This how- 
ever, does not happen because matching works so much 
better in ADMTT coordinates that the inner and near 
zone solutions are substantially closer. Also note that 
the singularities in the constraint violations of the data 
in ADMTT and harmonic gauge (see Figs. El and El 
occur at different coordinate locations. This is simply 
because the inner zone metric, which is relevant in this 
region, is expressed in different coordinates. 
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FIG. 16: Hamiltonian constraint violation of the global met- 
ric in ADMTT coordinates in the xy-plane for the equal mass 
system. The different line styles corresponds to isolines of 
constant constraint violation. 
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FIG. 17: This figure shows the y-component of the momen- 
tum constraint violation along the positive x-axis for the data 
computed in this paper (solid black line) and in paper 1 (dot- 
dashed black line) , assuming the parameters of the equal mass 
system. Note that all other components vanish along this axis. 



One obtains qualitatively similar plots if the con- 
straints are plotted along different directions. Evidence 
for this behavior can be seen in Fig. El which shows a 
contour plot of Hamiltonian constraint violation for the 
global metric in ADMTT coordinates in the xy-plane. 
We see that there is ring of radius 3m of negative con- 
straint violation around the black hole. Outside this ring, 
about 4.5m away from the hole, there are three maxima. 
The largest of these occurs on the x-axis. Note that these 
minima and maxima are all located in the buffer zone. In 
addition, there is a blow-up inside the horizon at ta = 0. 

In Fig. 1171 we plot the y-component of the momen- 
tum constraint along the positive x-axis close to black 
hole 1 for both data sets. For reasons of symmetry all 
other components vanish along this axis. Observe that 
the violation is everywhere small, reaching a maximum 
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FIG. 18: This figure shows the near zone (dashed line), the 
inner zone (dotted line) and the global lapse (solid line) along 
the positive x-axis in ADMTT coordinates for the equal mass 
system. Observe that the global ADMTT lapse crosses zero 
near the horizon. Also shown is the lapse (dash-dash-dotted 
line) in the new horizon-penetrating coordinates of Sec. IVIlH 



in the buffer zone. As expected, in the near zone the 
violation decays as we move away from the black holes 
and becomes identically post-Newtonian. Note however, 
that the resummed PN data in ADMTT gauge satisfy the 
momentum constraint exactly, while in harmonic gauge 
it has an error of 0(m^/r a) 2 '• On the other hand, in 
the inner zone, the violation is close to zero outside the 
horizon (which is located at ta ~ m^/2) and it grows 
as ta becomes larger. As in the case of the Hamilto- 
nian constraint the largest violation occurs in the region 
where the transition function leads to non-trivial aver- 
aging of the two approximations. The violation is once 
more smaller for the data presented in this paper relative 
to that of paper 1 , which is one more indication that the 
matching is smoother here. Finally, observe that close 
to the holes, and in particular close to the horizons, the 
violation due to the inner zone data diverges. This diver- 
gence can be traced back to the choice of slicing, which 
forces the lapse to be zero at the horizon. 

The vanishing of the lapse can be observed in Fig. 1181 
where we plot the inner zone, near zone and global lapse 
along the positive x-axis. This figure shows how the inner 
zone lapse goes through zero near the horizon. The van- 
ishing of the lapse translates into a divergence of the mo- 
mentum constraint, which renders excision difficult and 
the data difficult to implement numerically. However, 
this is a failure of the coordinate system used in the in- 
ner zone and not of the method of asymptotic matching. 
In the next section, we construct a new coordinate system 
in which the lapse remains non-zero across the horizon, 
and approaches an isotropic Schwarzschild lapse in the 
buffer zone to allow for matching. 

In Figs. El an d El we have compared the constraint 
violation of the data set presented here to that presented 
in paper 1. Technically this comparison suffers from the 
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problem that the constraints were computed in different 
coordinate systems. Note however, that both coordinate 
systems are identical to leading order. Yet instead of 
comparing the violations point by point, we can study 
their global average. Figs.ll4landll7lfavor the data of this 
paper because the humps are consistently smaller. As ex- 
plained earlier, this is an indication that after performing 
the matching, the inner and near zone solutions are very 
close to each other, making the transition smoother. 



VIII. HORIZON PENETRATING 
COORDINATES 

Recall that the global coordinate system constructed 
in Sec. [V] was based on matching the tidally perturbed 
black hole metric of the inner zone in isotropic coordi- 
nates to a post- Newtonian metric in ADMTT coordinates 
in the near zone. We have chosen to match in these co- 
ordinates because they are already very similar to each 
other. Hence, the coordinate transformation [Eq. (|26ll ] 
that leads to matching is close to the identity, and thus fa- 
cilitates computations. However, the fact that the global 
coordinate system remains very close to isotropic coordi- 
nates is also a disadvantage, since the t = slice is not 
horizon penetrating in isotropic coordinates. The lapse 
goes through zero very close to the horizon and in addi- 
tion, the extrinsic curvature has a coordinate singularity 
at the point where the lapse goes through zero. Initial 
data on such a slice may not be suitable for numerical 
evolutions. In this section we present a remedy for this 
problem. By constructing a further coordinate transfor- 
mation, we obtain coordinates which are close to horizon 
penetrating Kerr-Schild coordinates near each black hole. 

The basic strategy we will use is the following. We first 
determine the perturbed black hole metric valid in the in- 
ner zone in Kerr-Schild coordinates. We then transform 
this metric into a coordinate system which remains Kerr- 
Schild near the black hole horizons, but which is corotat- 



ing isotropic in the buffer zone. In this manner, the trans- 
formed inner zone metric components in the buffer zone 
will remain identical to those in Sec. El so that we can use 
the same matching coordinate transformation [Eq. il26[) ] 
as in Sec. [V] 

The standard transformation between spherical Kerr- 
Schild (t',f[) and spherical isotropic coordinates (i, f\) 
centered at hole 1 and its inverse are given by |44( 



t' = t + 2Miln 
i = i' - 2Mi In 



n 

2Mi 



V>(ri) 2 



2Mi 



V'(ri) 2 ' 



(47) 



where the radius is measured from the center of a black 
hole with mass M\ and where the conformal factor and 
its inverse are given by 



, / s M l 

<Kfi) = 1+2^, 



Mf') = 11. 



1-1- 



2Mi 



1/2 



(48) 



Note that the inverse transformation contains a square 
root in the conformal factor, which becomes complex in- 
side the horizon. For this reason, it is simpler to first 
transform the inner zone metric to Kerr-Schild coordi- 
nates analytically. Then, we can construct a new coor- 
dinate system that is the identity map near the horizon, 
but that brings back the metric to isotropic coordinates 
outside the horizon in the buffer zone. 

Let us first transform the inner zone metric of hole 
1 to Kerr-Schild coordinates. By applying Eqs. (|47|l to 
the inner zone metric in spherical isotropic coordinates 
[Eq. (3.22) of Ref. H3] we obtain 
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(49) 



where all coordinates are centered on BH 1 and where we 
used the abbreviation 

/ = (l - ■ (50) 

In Eq. (|49|) . i stands for the isotropic time coordinate, 
given in Eq. (|47|l . Recall that this metric was derived un- 
der the assumption that the second black hole responsible 
for the tidal perturbation is moving slowly. In particular 
Alvi obtained the perturbation from a stationary pertur- 
bation by replacing (f> by (f> — fit. This means the largest 
error in this perturbation is of order 0(m / b) 5 > 2 0(r a /b) 2 . 
In the following we will replace Schwarzschild time i by 
Kerr-Schild time t to simplify our expressions. This re- 



placement will change the tidal perturbation only at or- 
der 0(m/6) 5//2 0(r y 4/6) 2 and thus not introduce any extra 
errors. 

We now go one step further and transform Eq. Ij49(l to 
Cartesian Kerr-Schild coordinates via the standard map 

t' = t, f[ = (x 2 + y 2 + z 2 ) 1/2 , 

§ ' = C0S ^ / I 1/2 > 

(x 2 + y 2 + z 2 ) 
ft = tan" 1 (j^. (51) 

We then obtain 
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where ft = ft. This is the inner zone metric for hole 1 
in Cartesian Kerr-Schild coordinates, where the standard 
metric in Cartesian Kerr-Schild form is 



(KS) _ 



2M] 

ft 



(53) 



with null vectors l a — {l,x l /ft} and where we have in- 
troduced the shorthand 



3 [x cos fit + y sin fit) 



(54) 



The lapse in this coordinate system is now a positive 
definite function. 

The matching of Sec. E] however, is performed in the 
buffer zone in Cartesian corotating isotropic coordinates. 
We thus need a coordinate transformation that leaves 
Eq. I|52l) unchanged near the horizon, but takes the met- 
ric to isotropic corotating coordinates in the buffer zone. 
For black hole A, the transformation we use is given by 



t 


= (1 


x 


= (1 


y 


= (1 


z 


= (1 


where 






1pA 



= {l-Q A )t + Qj 



t + 2ttia In 



rA^A 
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1 



Qa) y + Qa [x sin fit + y cos 



fit) 7/4 



(55) 



rriA 
2f A ' 



f a = \/ 5? + y 2 + z 2 (56) 



and the transition function 

Qa = f(r A ) 

with 



To = 2.2777,4, 77J 

q = 1, s = 2.8, 



0.6r¥ - 2.2m A , 



(57) 



(58) 



is designed such that Qa is unity in the buffer zone and 
zero near black hole A (see Fig.[5J|. This means that using 
Eq. (|55|l we can transform from Kerr-Schild coordinates 
(labeled by a hat), to coordinates (labeled by a tilde) 
which are corotating isotropic coordinates in the buffer 
zone (where Qa = 1) ; but which are equal to Kerr-Schild 
coordinates at the black hole horizons (where Qa = 0). 
The function Qa is chosen carefully so that the trans- 
formed coordinate system is identically Kerr-Schild at 
the horizon (ta — 2mA), while it is equal to isotropic co- 
ordinates in the buffer zone where we perform the match- 
ing. We can adjust, as usual, how fast we transition by 
changing the parameters of Qa, but we are constrained 
to having the metric completely in isotropic coordinates 
in the buffer zone, in order for the matching to be valid. 
The inner zone metric then becomes 



o (1) J?J /3 



(59) 



where, in analogy to Eq. I|31|) . the Jacobian matrix is 
given by 



J H 



dx & 
dxi 1 



(60) 
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FIG. 19: Schematic drawing of the different metrics, coor- 
dinate systems and transition windows. The solid line rep- 
resents the distance r\ (in ADMTT coordinates) from BH 1. 
The thick solid line labeled by ADMTT delineates the region 
where we use the PN near zone solution in ADMTT coor- 
dinates. The other two thick solid lines labeled by ISO and 
KS, show the region in which we use the inner zone solution. 
The labels here stand for the isotropic (ISO) and Kerr-Schild 
(KS) coordinates we have used in the inner zone before apply- 
ing the matching coordinate transformation 12611 . The dotted 
line Tl corresponds to the region where the transition func- 
tion Qi (see Fig.|SJ of Eqs. 15711 and 1581 is between zero and 
one, so that the coordinate transformation of Eq. 155H yields 
non-trivial results. I.e. in that region the inner zone metric 
is given in coordinates that are a mixture of Kerr-Schild and 
isotropic coordinates, while to the left and right of Tl it is 
given in Kerr-Schild and isotropic coordinates, respectively. 
Using the definitions 15711 and 1581 1 for the transition function 
Qi , we find that the inner and outer radii of Tl are 2. 2m a and 
0.6r^, where r A ' given in Eq. The dotted line labeled 

T2 indicates the transition window of the transition function 
Fi (see Fig.|5J of Eq. 14211 with parameters given by Eqs. 14311 
and 139H . In T2 we use a weighted average of both the in- 
ner zone metric (transformed into ADMTT coordinates using 
Eq. lH3) and the near zone ADMTT metric. From Eq. 
we see that the inner and outer radii of T2 are 0.4rjf and 
3.9r¥. 



This Jacobian can be computed by taking derivatives of 
Eq. JS3>- 

With this new coordinate transformation, the global 
metric contains several transition functions. A schematic 
drawing of these transitions is presented in Fig. 1191 In 
this drawing, the black dot represents BH 1 and the solid 
black line is the radial direction, where the horizon, for 
example, is located at f\ = 2m\. Fig. ^| also shows the 
different coordinate systems used, where KS stands for 
Kerr-Schild, ISO for isotropic coordinates and ADMTT 
for the PN near zone. The dotted line shows the region 
where the transitions take place: Tl is the transition be- 
tween KS to ISO and T2 is the transition produced by 
the matching coordinate transformation. Observe that 
Tl is chosen so that the 4-metric is in KS coordinates 
everywhere near and inside the horizon, while it is com- 
pletely in ISO coordinates where T2 begins. This restric- 
tion makes the transition window of Tl narrow. On the 
other hand, the transition window of T2 is restricted only 
by the size of the buffer zone and, thus, is chosen to be 
wider. 



Effect of the new coordinates 

In this subsection, we describe the advantages of the 
new coordinate system. In these coordinates, the new 
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FIG. 20: This figure shows the xx-component of the global 
metric in new coordinates for an equal-mass system along the 
x-axis. 



global 4-metric is much better behaved close to the hori- 
zons, where it is of Kerr-Schild form. Thus, the lapse re- 
mains positive definite through the horizon. The lapse for 
an equal mass system along the positive x-axis in the new 
coordinates is shown in Fig. ^| (dash-dash-dotted line) 
together with the lapse in the old ADMTT coordinates 
(solid line). The inner zone lapse in new coordinates 
is approximately equal to the Kerr-Schild lapse inside 
the horizon, while it smoothly approaches the lapse of 
isotropic coordinates in the buffer zone and the ADMTT 
lapse in the near zone. 

Another quantity that changes in the new coordinates 
is the spatial metric. Fig. 121)1 shows the ^-component of 
the new global metric along the x-axis for an equal mass 
system. The large humps on either side of each black 
hole in Fig.|2U|are produced by derivatives of the coordi- 
nate transformation (|55|l . which contains the transition 
function Qa that changes rapidly from zero to unity in 
a small region of width w = O.Qr^ — 2.2mA- Had we 
chosen a wider transition window w to transition, this 
derivative would become smaller and the hump would be 
smoothed out. However, we need the metric to be com- 
pletely in isotropic coordinates in the buffer zone where 
we perform the matching (approximately at ta — Qhtia)- 
Therefore, we are constrained to have a narrow transi- 
tion window w, which then produces large derivatives of 
the transition functions and humps in the metric. Note, 
however, that these humps are not spurious gravitational 
radiation. They simply arise because of performing a 
coordinate transformation. Therefore, due to the inher- 
ent diffeomorphism invariance of General Relativity, the 
physical content of the data will not be altered by the 
coordinate transformation. In addition, if we choose a 
larger black hole separation, increases and we obtain 
a wider transition window to transition, and hence the 
humps become smaller. 

Since the inner zone lapse in this new coordinates is 
now a positive definite function, we expect the extrinsic 
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FIG. 21: Plot of the y-component of the momentum con- 
straint violation for an equal mass binary near BH 1 in new co- 
ordinates (black solid line) and in ADMTT coordinates (dot- 
dashed line). Observe that in new coordinates there is no 
divergence until the singularity. 



FIG. 22: This figure shows the Ricci scalar at t = along 
the :r-axis for the inner (dotted line), near (dashed line) and 
global (solid line) 4-metrics. 



IX. CONCLUSIONS 



curvature and the momentum constraint violation to be 
small and well-behaved across the horizon. In Fig. 1211 we 
compare the y-component of the momentum constraint 
violation for an equal mass binary in this new coordinate 
system to that in the old ADMTT coordinates near BH 
1. Observe that the constraint violation is finite every- 
where, except near the singularity. With this new coordi- 
nate system, excision is now possible, since the curvature 
does not blow up until close to the physical singularity. 
We have zoomed to a region 0.5m away from the singu- 
larity of BH 1 to distinguish the behavior of the violation 
better. As we can observe from the figure, the ADMTT 
constraint violation is identical to the violation in the 
new coordinates away from the horizon. However, near 
the horizon there are spikes when we use the old coor- 
dinates. These spikes are poorly resolved in this figure, 
but we have checked that they are indeed divergences. 
Observe that these spikes are not present in the new co- 
ordinates. 

Since our method really yields a 4-metric we can also 
compute the Ricci tensor. This tensor should be close 
to zero for an approximate vacuum metric. Figure 1221 
shows the Ricci scalar along the x-axis at t = 0. One 
can see that it has the same qualitative features as the 
Hamiltonian and momentum constraints. I.e. apart from 
the singularity at the center of the black hole, the largest 
violations occur in the buffer zone. Notice however, that 
unlike in the case of metric (see Fig. I20|l the coordi- 
nate transformation 1)55(1 ■ with the transition function 
Qa does not produce large humps in the Ricci scalar. 
This confirms that the large humps in the metric are a 
pure coordinate effect, that has no influence on coordi- 
nate invariant quantities such as the Ricci scalar. 



We have constructed approximate initial data for non- 
spinning black hole binaries by asymptotically match- 
ing the 4-metrics of two tidally perturbed Schwarzschild 
solutions in isotropic coordinates to a post-Newtonian 
metric in ADMTT coordinates. The two perturbed 
Schwarzschild metrics are valid only close to each black 
hole (in the inner zones) and contain tidal deformations 
which are correct up to 0(r A /b) 2 . The ADMTT post- 
Newtonian metric we use is formally correct only up to 
0(m/6) 3//2 and valid only in the near zone, and not close 
to the black holes. However, by adding certain higher 
order terms to the ADMTT metric we bring it into a 
form that is close to the Schwarzschild metric plus some 
artificial perturbation. This means that instead of blow- 
ing up in an uncontrolled way, the ADMTT metric ap- 
proaches the Schwarzschild metric in isotropic coordi- 
nates when we approach one of the black holes. Since 
the two tidally perturbed Schwarzschild solutions, which 
are actually valid close to each black hole are also given 
in isotropic coordinates, both metrics agree at leading 
order even before asymptotic matching is performed. 

The procedure we have used to achieve asymptotic 
matching closely follows the calculation of paper 1, and 
is performed in the so called buffer zone where both 
the post-Newtonian and the perturbed Schwarzschild ap- 
proximations hold. The result is that both metrics agree 
in the buffer zone, up to the errors in the approxima- 
tions. However, since both approximations are similar to 
Schwarzschild in isotropic coordinates, matching yields 
much better results than in paper 1, where harmonic co- 
ordinates were used for the post-Newtonian 4-metric. In 
addition, after matching in the buffer zone the two met- 
rics are also very similar near the black hole horizons, 
even though the post-Newtonian metric is formally not 
valid there. The biggest deviation of the post-Newtonian 



24 



4-metric from the correct perturbed Schwarzschild solu- 
tion near the horizon occurs in the <?oi or shift compo- 
nents at 0(m/6) 3 / 2 . At this same order similar devia- 
tions are also observable in the extrinsic curvature of the 
t = slice. 

The resulting global piece-wise 4-metric is made for- 
mally C°° by the use of smoothing functions. These 
functions are chosen such that the smoothed global met- 
ric has errors on the order of the error introduced by the 
more accurate of the two approximations we match. This 
smoothing procedure works much better than in paper 
1, because the tidally perturbed Schwarzschild solutions 
and the post-Newtonian metric are much closer to each 
other and diverge at similar locations. The smoothed 
global metric is obtained in ADMTT coordinates, and 
thus very similar to isotropic coordinates near each black 
hole. Hence our coordinates are not horizon penetrat- 
ing and, for example, result in a lapse that goes through 
zero close to the horizons. Since such coordinates may 
be problematic for certain numerical simulations, we con- 
struct an additional coordinate map. With this map one 
can transform the 4-metric from global ADMTT coor- 
dinates, obtained through asymptotic matching, to new 
coordinates which are similar to Kerr-Schild coordinates 
near each black hole, but which remain ADMTT further 
away from the black holes. These new coordinates are 
horizon penetrating and lead to a lapse which is every- 
where positive on the t' = slice. The same map can in 
principle also be applied to the data computed in paper 
1. 

For both, global ADMTT and new coordinates we have 
constructed initial data sets on the respective t' = 
slices. These initial data are then used to compute Hamil- 
tonian and momentum constraint violations. We find 
that the initial data found in this paper is closer to the 
constraint hypersurface that the data presented in pa- 
per 1, making it perhaps more amenable to numerical 
simulations. 

In conclusion, we have used the method of asymptotic 



matching to construct improved approximate initial data 
for non-spinning black hole binaries. Future work will 
concentrate on applying this method to spinning binary 
systems. This problem is of fundamental importance to 
astrophysics and gravitational wave detection, since most 
black holes are believed to be rotating. Using the results 
of the matching of a tidally perturbed Kerr black hole 
to a post-Newtonian metric should be possible. Other 
work will concentrate on repeating this analysis to higher 
post-Newtonian order to explicitly incorporate the effect 
of gravitational waves. 
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APPENDIX A: JACOBIAN 

In this section we present explicit formulas for the Ja- 
cobian of the transformations given in Eq. (|26l) . These 
formulas are long and complicated, but, if the reader is 
interested in implementing them, we can provide a Maple 
file with explicit expressions for them. The Jacobian is 
given by 
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where terms not listed here are zero. 

This is the Jacobians of the transformations that allows 
us to construct a global metric by transforming the inner 
zone metrics in the buffer zone. Note that this is the 



Jacobian for the transformation that is valid in buffer 
zone 1 (0i3.) In order to obtain the Jacobian for the 
transformation in the other buffer zone (O23), one can 
apply the substitutions of Eq. (|2*5|l . 
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